Random Exchange Quantum Heisenberg Chains 



in 



A. FurusakiJ] M. Sigrist, E. Westerberg, P. A. Lee 
Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139 

K. B. TanakaB and N. Nagaosa 
Department of Applied Physics, University of Tokyo, Hongo, Bunkyo-ku, Tokyo 113, Japan 

(February 1, 2008) 

The one-dimensional quantum Heisenberg model with random ± J bonds is studied for S — \ and 
S = 1. The specific heat and the zero-field susceptibility are calculated by using high-temperature 
series expansions and quantum transfer matrix method. The susceptibility shows a Curie-like tem- 
perature dependence at low temperatures as well as at high temperatures. The numerical results 
for the specific heat suggest that there are anomalously many low-lying excitations. The qualitative 
nature of these excitations is discussed based on the exact diagonalization of finite size systems. 
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I. INTRODUCTION 

One-dimensional (ID) quantum spin systems have at- 
tracted much interest for many years. They provide a 
very useful test ground for our understanding of quantum 
systems with many degrees of freedom and have led to 
the development of various powerful methods of theoreti- 
cal physics. Quantum spin systems are also of interest in 
their own right for the large variety of their ground states 
and properties of excitations. Over the years the search 
for real spin chains has led to the synthesis of various 
quasi-lD systems. In general we encounter defects and 
disorder in real systems, which may destroy much of the 
properties expected for uniform spin chains. Although at 
first sight unwanted, these defects have initiated a great 
deal of studies on disordered chains of classical as well as 
quantum spins. 

To our knowledge the first system investigated 
in this respect originates from the quasi-lD or- 
ganic charge transfer compounds of tetracyanoquin- 
odimethanide (TCNQ) which was described by the 
nearest-neighbor Heisenberg model, 
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where the interaction is antiferromagnetic but random 
in strength with a continuous broad probability distri- 
bution, P(Ji)B This type of model was investigated by 
Hirsch and Josecl by means of renormalization group 
(RG) methods based on a Kadanoff block spin decimen- 
tation scheme. Ma and coworkers,!] used a different real 
space RG method where the strongest bonds were suc- 
cessively integrated out in order to reach an effective 
Hamiltpnian with a fixed point of the distribution P(Ji). 
Fishery recently revisited this RG method and extended 
it also to the XXZ Heisenberg model. The result of 
this analysis is that the ground state is a random sin- 
glet phaseO where each spin is paired into a singlet with 
another spin which may be located far away in the chain. 



Another approach for the disordered spin chains is the 
mapping of the spins into fermions via a Jordan- Wigner 
transformation, which in turn are treated by a bosoniza- 
tion scheme. This approach was used by NagaosaQ to 
study the effect of a random magnetic field. He de- 
rived scaling relations for spin-spin correlation functions 
analytically and confirmed them by numerical calcula- 
tions. Assuming weak disorder, Doty and Fisheru then 
analyzed the phase diagram and thermodynamic prop- 
erties of XXZ spin chains with various type of disorder 
by applying a perturbative RG method to the bosonizcd 
Hamiltonian. These results have been confirmed to some 
extent by extensivej-pumerical calculations for finite sys- 
tems by Haas et alB and by Runge and Zimanyia 

A new class of random spin chains has been found in 
the compound SraCuPtOe where the Cu-site alternates 
with the Pt-site along chains .1121 The Cu has a spin-i de- 
gree of freedom while Pt has none. The interaction be- 
tween the Cu spins is antiferromagnetic. Measurments 
of the uniform susceptibility show that this system is 
well described by a ID antiferromagnetic (AF) Heisen- 
berg model. The Pt atoms can be replaced by Ir, which 
carries a spin i like Cu. In the chains the spin on an 
Ir atom couples ferromagnetically to its neighboring Cu 
spins. As a result the number of spins has doubled and 
the interaction is ferromagnetic. This spin chain is de- 
scribed well by the ID ferromagnetic (FM) Heisenberg 
model. It is also possible to produce an alloy by replac- 
ing a fraction 1 — x of P|t-by Ir: SraCuPt^Iri-^Og, where 
x can be chosen at will Jill This yields a system with ran- 
dom FM and AF interaction of fixed strength. 

Motivated by these exp eriments we consider the 
Heisenberg model, Eq. (1.1), where the probability dis- 



tribution of Ji has the following form, 

P{.h) = P 5(Ji + J) + (1 - p)6(J l - J). 



(1.2) 



The coupling J^'s have the same strength, but are ran- 
dom in sign, with a probability of p (0 < p < 1) to be 
ferromagnetic and f — p to be antiferromagnetic. We 
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note that this model does not exactly realize the con- 
ditions found in the alloys described above. In the real 
system the strengths of the FM and the AF bonds are 
in general different. Additionally, there is the correlation 
(neglected in our model) that FM bonds appear always 
in even numbers. Nevertheless, we believe that the char- 
acteristics of our model are very similar to the ones of 
the real alloy and that the properties which determines 
the physics is the discrete randomness with FM and AF 
bonds. 

Obviously the analytical treatments mentioned above 
are not immediately useful for our model. The scheme 
by Ma and coworkers cannot be simply applied here, be- 
cause there are no strongest bonds which would allow the 
decimation procedure. Furthermore, the disorder is not 
weak in our case, and thus perturbative treatments would 
not lead to sensible results. Instead we apply various 
types of numerical methods, such as high-temperature 
expansion (HTE), transfer matrix method (TMM), and 
exact diagonalization (ED) of small systems. The first 
two provide accurate results over a wide range of tem- 
perature (k B T > 0.1J) and the last one can be used to 
reach a qualitative understanding of certain aspects of 
the low-temperature regime. 

Before going into technical details we will briefly out- 
line the picture which emerges from our analysis. The 
main result is that this quantum spin system exhibits 
three distinct regimes, high-, intermediate- and low- 
temperature regime. It will become clear that this prop- 
erty originates from the fact that the spin chain consists 
of alternating FM and AF segments of variable length. 
At very high temperatures (&bT ^> J) the interaction 
among the spins does not play any role and the spins 
behave independently. Hence, the randomness of the 
exchange does not play any role in the physical prop- 
erties. In particular, the uniform susceptibility follows 
the ordinary Curie law, \ 00 When the tempera- 

ture is lowered, the spins gradually start to correlate. In 
the AF segments the spins form collectively the smallest 
possible total spin (either a singlet, S = 0, or a dou- 
blet, S — |), while the FM segments build up a large 
spin degree of freedom by aligning spins parallel to each 
other. Due to the quantum nature of the spins the in- 
teraction among these new effective degrees of freedom is 
rather weak. Therefore, in an intermediate temperature 
range (k B T < J), they behave like independent spins 
and lead to a new Curie behavior of the susceptibility, 
however with a different (effective) Curie constant. The 
excitations relevant in this temperature regime are dom- 
inated by the discrete spectrum of "spin wave" modes 
of each segment. These intrasegment modes are local- 
ized because of the mismatch of momentum and energy 
between the modes in different segments. Thus these 
modes cannot be transfered easily between different seg- 
ments. The thermodynamics of this regime is effectively 
described by an ensemble of decoupled FM and AF seg- 
ments of finite length. Only at very low temperatures 
(fcflT <C J) the interactions between spins belonging to 



different segments become important. Within our treat- 
ments this low-temperature regime is least understood 
at present. However, our results suggest that a consider- 
able fraction of density of states is located at low ener- 
gies (u -C J). It is characterized by an effective model 
of spins with various sizes which are coupled by cither 
ferromagnetic or antifcrromagnetic interaction of widely 
distributed strengths. 

The crossover between the regimes manifests itself in 
peak-like structures of the specific heat. These structures 
are located where spins begin to correlate. While the uni- 
form FM and AF systems have only one such (broad) 
peak (corresponding to the correlation of the S = | 
spins), we expect that two structures appear in the ran- 
dom system, a peak at temperature where the original 
S = | spins start to correlate (k B T ~ J) within each 
segment, and, possibly, a shoulder where the effective 
spin degrees of freedom in each segment start to corre- 
late. The location and size of the second structure will 
depend on p. 

The existence of three regimes is in contrast to the 
situation in the analogous classical case, where only two 
regimes exist. The reason for this difference lies in the fol- 
lowing properties of the classical spin system. In the case 
of classical spins the bond disorder of the type of Eq. ( |l.2j ) 
is irrelevant for the thermodynamic properties. A simple 
transformation of the spin variables (changing the sign 
of some of the spins) yields a uniform FM spin chain 
without affecting the energy spectrum. Consequently, 
the specific heat would not depend on the disorder and 
occurrence of FM and AF segments does not have any 
direct implications here. However, bond disorder affects 
the spin-spin correlation and the susceptibility. The sus- 
ceptibility x shows a clear signature of a crossover from 
a high-temperature regime to a new regime at lower tem- 
perature similar to that seen for the quantum spin sys- 
tem. Despite this qualitative similarity the nature of the 
crossover is quite different. The classical spins do not sep- 
arate into segments, but they correlate within a length 
£, the correlation length of the uniform system, which is 
independent of disorder. As this length grows with de- 
creasing temperature, the spins in a cluster of length £ 
act together as one large effective spin whose size scales 
as S e g oc i; 1 / 2 for £ much larger than the lattice constant. 
(This follows from the fact that the directions of the spins 
are random depending on the signs of the bonds so that 
a random walk picture applies). These effective spins be- 
have independently and the susceptibility per lattice site 
shows the Curie behavior in the low-temperature limit: 



It is obvious that there are no further correlation ef- 
fects beyond this, so that this second Curie regime repre- 
sents indeed the low-temperature regime reaching down 
to T = OK. 

This paper is organized as follows. In Sec. II we con- 
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sider technical aspects of the HTE and the TMM, in 
particular, the problem of extrapolation from the high- 
temperature limit to low temperatures. Readers who are 
not interested in technical details can skip this section. 
In Sec. Ill we present the numerical results of the HTE 
and the TMM. Section IV is devoted to the results of 
the ED of small clusters. We demonstrate the weakness 
of the coupling between the effective spins formed in the 
segments. Finally, we discuss our results in Sec. V. 



II. METHODS OF NUMERICAL ANALYSIS 

In this section we describe some technical details of the 
HTE and the TMM. The former allows for an exact en- 
semble average over the random exchange for arbitrary p. 
We use the latter method only for p = 0.5 to check the 
validity of the extrapolation scheme used in the HTE. 
It turns out that both methods give good quantitative 
results over a wide temperature range kgT > 0.1 J. 



A. High-temperature series 



We then average over all the configurations for given 
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In the next step we recursively subtract the contributions 
of all subclusters from o n (K;p): 

n-l 

o' n {K-p) = o n (K;p) -J2(n-k+ l)o' k (K;p). (2.4) 
fc=i 

Note that this subtraction is only possible after the av- 
erage has been taken. Summing up the series o' n (K;p) 
yields the final series 



o(K-p) 



N 



o' n (K;p) 

^^ G,m)^y™. 



(2.5) 



The high-temperature expansion of spin systems can 
be implemented as a-jcluster expansion algorithm in a 
straightforward way.Ea In ID this algorithm is greatly 
simplified. 

Consider a cluster of length n which is described by 
Hamiltonian, 
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(2.1) 



and assume that m of (n — 1) exchange couplings, Ji, 
are ferromagnetic and the others are antiferromagnetic. 
There are n C m = n\/m\(n — m)\ different configurations 
of this type. For each configuration of couplings {Ji} 
we calculate thermal average of a physical quantity in 
powers of the inverse temperature (3: 
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(2.2) 



where K — ^f3JS. The thermodynamic quantities we 
consider are the internal energy u and the uniform sus- 
ceptibility x, for which the operator O n is O n = TL n and 
On = /VQXi 5 ?) 2 <> : Bohr magneton). The coeffi- 
cients di iTntn ({ Ji}) have to be evaluated for all configura- 
tions {Ji}. However, some operations (such as reflection 
in the center of the cluster and change of the sign of Ji) 
allow us to reduce the number of configurations we need 
to calculate. 



For ID systems it t urns out that the high-temperature 
series given by (2.5) is correct-up to the order K 2N for 
u, but only up to K N for x.EJ The largest system size 
we considered is N = 11 (N — 8) for the spin-i (1) case. 
In Tables §-|y] we give lists of the expansion coefficients 
u(l, m) and x(h m )- For the spin-i AF chain (p = 0) o 
result agrees with the earlier calculation by Baker et al 



B. Analysis of the series 



The high-temperature series of finite length by itself 
does not give a reliable result for ksT < J. There 
are, however, various extrapolation schemes which can 
provide a very—good approximation down to rather low 
temperatures.E3 Since no finite-temperature transition is 
possible for our ID spin systems, one should, in principle, 
be able to extrapolate the high-temperature series down 
to zero temperature. 

For our analysis we use the Pade approximation 
method.E.3 Although at first sight it seems straightfor- 
ward to apply it to the series of the internal energy 
u(K;p), we encounter the following severe problem, 
which comes from the very fact that the series could be 
extrapolated down to zero temperature or K — > oo. From 
the relation u(K;p) = —u(—K; 1 — p), we see that the 
series u(K; x) for the spin chain with p = x also describes 
the spin chain with p = 1 — x on the negative real axis. 
In particular, u(— oo;p) is a modulus of the ground state 
energy of the latter system while u(oo;p) is the ground 
state energy of the former. That is, the series must have 
different limits for K — » ±oo. It is impossible to deal with 
this feature in the standard Pade approximation scheme. 
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To circumvent this difficulty, we use the following 
method. We first note that u(K;p) is a monotonic func- 
tion of K along the real axis. Thus we can invert the 
series u(K;p) to K(u;p), a series in powers of u/ J S3 The 
function K(u;p) vanishes at u — and diverges both at 
u = u(oo;p) and at u = — u(oo; 1 — p). Thus the tem- 
perature range — oo < K < oo is mapped to the region 
u(oo,p) < u < —u(—oo; 1 — p). 

The inverted series K(u;p\-can now be analyzed by 
means of Pade approximantstj because the two limits 
K — > ±oo correspond to two finite points. If we assume 
that near zero temperature the internal energy behaves 
as u(K;p) — u(oo;p) oc T 1+Q , or equivalently 



K oc [u — u(oo,p)] 



(2.6) 



then the ground state energy u(oo;p) should show up in 
^-/K as a pole on the real axis. Thus we can estimate 
the ground state energy and the exponent a from Pade 
approximants of ^-/K. The ground state energy we es- 
timated from these Dlog Pade approximants is given in 
Table We expect that the estimated numbers have 
errors of order 1% except p = 1, for which the Pade ap- 
proximants have disturbing poles close to the physical 
one corresponding to the ground state energy. 

It turns out, however, that we cannot obtain reliable 
estimates for the exponent a in this way; the estimated 
values seem to be too large. For example, for p = Pade 
approximants give a ~ 1.85, which differs significantly 
from the exact value a = 1. This failure in estimating 
a is possibly related to the fact that the singularities 
at it = u(oo;p) and — u(oo; 1 — p) in the Dlog Pade ap- 
proximants are not simple poles. In fact, the approxi- 
mants have many poles on the real axis, suggesting that 
there are two cuts along the real axis: u < u(oo;p) and 
u > — u(oo; 1 —p). We suspect that the exponent is more 
sensitive to the disturbing cuts than the location of poles 
(ground state energy) is. 

Now we are ready to explain how the specific heat 
and the susceptibility are calculated from the high- 
temperature series. By integrating the Pade approxi- 
mants for ^-/K, we first obtain K(u;p). The specific 
heat is then given by 



C{u(T))=-[K{u)\- 



du 
dK' 



(2.7) 



Due to the too large a, the specific heat calculated from 
Eq. ( 2/7 ) is too small at very low temperatures. As we 
will see in the next section, however, it still gives reason- 
able results for fegT > O.IJ. 

The susceptibility is calculated from the high- 
temperature series x(K]p) m the following way. We first 
transform the series of Tx, which is a power series of K, 
to a new series in powers of u by substituting the inverted 
series K(u) into T\{K ). We then make a Pade approxi- 
mant for -^Tx{u) /Tx(u) and integrate it back to get an 
approximant for Tx(u), which we denote P(u). The sus- 
ceptibility is finally calculated from x(T) = (3P(u(T)). 



The results obtained in this way are shown in the next 
section. 



C. Quantum transfer matrix method 

In addition to the high-temperature series expansion, 
we use the quantum transfer matrix method, which is 
also a numerical method valid at high temperatures, to 
check the validity of the extrapolation scheme used in 
analyzing the high-temperature series. The basic idea of 
the method is to calculate the partition function by using 
the Trotter breakup, 

exp(-/3W) w [exp(-/3H c /iVT)exp(-/3H /iVT)] ArT , 

(2.8) 

where TL = Ji c + H. with a checkerboard decomposition. 
Then the system can be viewed as a two-dimensional clas- 
sical system, and the partition function is calculated by 
multiplying a transfer matrix. The details of the method 
can be found in the liter atures.EJ This method has the 
following advantages: 1) The CPU time grows only lin- 
early with the system size L so that we can study large 
systems. The largest system we studied has 400 sites, 
for which we can expect self-averaging. 2) Although the 
Trotter number is limited by the memory size of the 
computer [Nt — 10 (6) for S — h (1) in our calculation], 
it is possible to extrapolate to the Nt — > oo limit—from 
the free energy Fn t calculated for a finite Nt, byE3 



ai 
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N T 



(2.9) 



For example, for the spin-i case we use Fiq, Fg, and 
Fg to determine F^, a±, and a,2- This method works 
well down to kgT ~ O.IJ. 3) In contrast to the quan- 
tum Monte Carlo method, there is no statistical error in 
the transfer matrix method. Thus, one can calculate the 
specific heat and spin susceptibility by numerical differ- 
entiation with respect to the temperature and a magnetic 
field. 



III. THERMODYNAMIC PROPERTIES 

In this section we discuss the results for two measur- 
able quantities, the specific heat and the susceptibility, 
for the spin-i and spin-1 chains. 



A. Specific heat 

As described in Sec. IIB, the specific heat is calcu- 
lated from the high-temperature series by Eq. (2.7). We 
show the results for p = 0, 0.25, 0.5, 0.75, and 1 in Figs. 

Note that the results for the spin-i AF chain 



|T| and 
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(p = 0) and for the FM chain (p — 1) are due to two-point 
Pade approximants which we impose to have poles at the 
known ground state energy with the correct exponent a 
(a = 1 for p = and a = 0.5 for p = 1). The quality 
of our approximants is demonstrated by the comparison 
with .the data obtained from finite size calculations by 
Biotech for p = 0, 1 and by the data we obtained using 
the TMM for p — 0.5. We see that our Pade approxi- 
mants are quite reliable at least above k B T ~ 0.1J for 
5 = | and k B T ~ 0.2 J for S = 1. 

However, we find that our approximants for < p < 1 
cannot be valid in the whole temperature range; they do 
not satisfy the sum rule, 



fi 2 S 2 
Sk B T 



Sac — 



C(T) 
T 



dT = k B ln(25 + l). 



(3.1) 



This can be clearly seen from Figs. |^ and ^, where we 
show C(T)/T as a function of T. Some entropy is miss- 
ing if we naively extrapolate our high-temperature se- 
ries down to zero temperature. Rough estimates give 
the following values for the missing entropy: in the spin- 
et (spin-1) system AS/S^ =12%, 24%, and 36% (14%, 
28%, and 32%) for p =0.25, 0.5, and 0.75, respectively. 
Since our results are reliable above k B T ~ (0.1J — 0.2J), 
we expect that-the missing entropy is "hidden" at lower 
temperaturesO'Ej A considerable fraction of the density 
of states would be located at low energies. We discuss in 
Sec. IV the nature of the low-lying excitations which are 
responsible for the missing entropy. 



B. Susceptibility 

We plot the susceptibility calculated from the HTE in 
Figs. H and [| For p = 0.5 we also plot the data obtained 
from the TMM. The results obtained from the two dif- 
ferent methods agree quite well so that we can trust our 
data down to low temperature ~ 0.1 J. We note that the 
convergence of the data of the susceptibility in the TMM 
is much better than that of the specific heat. 

From Figs. || and ^ we see that the susceptibility obeys 
the Curie law at low temperatures as well as at high tem- 
peratures. There is a crossover at k B T ~ J, and the 
Curie constant at lower temperature depends on p. In- 
terestingly enough, this feature is qualitatively the same 
as that of the classical Heisenberg spin chain: 



Xd(T) 



li 2 S 2 l + v{JS 2 /k B T) 
3k B T l-v{JS 2 /k B T) 



with 



v(x) = (2p- 1) 



coth(2a;) — — 



(3.2) 



(3.3) 



Equation— (3.2) is obtaineaSEH by using Fisher's 
method.E3 For < p < 1 the susceptibility shows a Curie- 
like behavior for both T — > oo and T — > 0: 



Xd(T) 



for T 



and 



Xd(T) = 



p 2 S 2 p 
3k B T 1 — p 



for T -> 0. 



(3.4) 



(3.5) 



Obviously there is a crossover between the two regimes at 
k B T ~ J where the Curie constant changes its high- to 
its low-temperature value. Note that at p = 0.5 the sus- 
ceptibility of the spin chain is the same as that of a free 
spin. It is also worth noting that for the uniform system 
the low-temperature susceptibility is not Curie-like, but 
goes to a constant for p = and diverges quadratically 
for p = 1, which is qualitatively the same behavior as the 
spin-i chains. Ej 

Besides the qualitative similarities, there are clear dif- 
ferences between the quantum and the classical spin 
chains. As discussed in Sec. I, the qualitative nature 
of the crossover is completely different. For the classical 
spin chains the existence of the length scale £ is essential, 
whereas the segmantation is crucial for the quantum spin 
chains, as will be demonstrated in the next subsection. 
A quantum effect is clearly visible for p = 0.5 where, 
in contrast to the classical case, we find a crossover at 
k B T ~ J between the high-temperature Curie law, 



X 



li 2 S{S + l) 
3k B T 



(3.6) 



and a low-temperature Curie-like behavior. For larger 
spins, however, the change becomes less pronounced as 
we can see in Fig. ^. 

Unfortunately, we cannot simply claim from our nu- 
merical result that the susceptibility diverges as 1/T for 
T — > like in the classical system. Although we are con- 
fident about the accuracy of our data above k B T ~ 0.1J 
from the comparison between the HTE and the TMM, 
we cannot assume that our extrapolation scheme is valid 
down to zero temperature, which is certainly not the case 
for the specific heat. To determine the T-dependence of 
the susceptibility at T — > 0, we need to develop a com- 
pletely different approach which can deal with low-energy 
excitations directly. 



C. Discussion of the results 

In the introduction we mentioned that the key to the 
understanding of our random quantum spin system lies 
in the property that the spin chain consists of a sequence 
of alternating AF and FM segments of various lengths. 
On an intermediate energy scale they behave like inde- 
pendent finite size systems. Thus, for low enough tem- 
perature the spins within each segment correlate into the 
ground state as if the segments were decoupled. There 
are small thermal fluctuations to the excited states which 
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are separated by a finite energy gap because of the finite 
length of each segment. The ground state of the AF seg- 
ments is a spin singlet or doublet (triplet) for an odd or 
an even number of bonds in the spin-i (spin-1) system. 
In the FM segment the ground state is characterized by 
the formation of the largest possible total spin, where 
the boundary spins should not be included, i.e., the total 
spin of a segment with £ bonds is Stot = S(£ — 1). It is 
important to realize that the boundary spins "shared" by 
adjacent FM and AF segments do not contribute to Stot 
of the FM segment, but always have to be counted to the 
AF segment for the formation of its ground state. This 
fact is confirmed by analytic arguments and finite size 
calculation, which will be dicussed in the next section. 

For short segments the intrasegment excitation ener- 
gies are of the order of J. For segments of larger length 
£ the lowest excitation energy can be estimated based on 
a spin wave picture in a finite chain, which scales as 



AE 



.]_ 

J 
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in FM segments, 
in AF segments. 



(3.7) 



The intrasegment spin-wave-like modes are localized 
within each segment, because the excitation spectra in 
neighboring segments are different. Therefore we expect 
that also the coupling between the effective spin degrees 
of freedom (formed in the ground state of each segment) 
is very weak. On an intermediate energy scale AE) 
they may be considered as independent spins. Only for 
much lower energies correlation would develop among 
them as will be demonstrated in the next section. 

These properties allow us to interpret some of our 
results of the HTE and the TMM. Let us first con- 
sider the specific heat. As shown in Figs. |l| and [|, 
the specific heat of disordered spin chains (0 < p < 1) 
drops more rapidly at low temperatures than that of uni- 
form chains (p = 0, 1). As a consequence the entropy 
we obtained numerically does not satisfy the sum rule, 
Soo = fcs ln(2S + 1) per site. From our discussion it be- 
comes clear that the effective spin degrees of freedom do 
not correlate in the temperature range covered by the 
HTE. Therefore their contribution to the entropy is not 
visible in the extrapolated specific heat. Instead they 
would appear at lower temperature which is beyond reach 
of our methods. A simple estimate of the missing entropy 
can be given by regarding all the segements to be inde- 
pendent. In this picture a FM segment contributes the 
entropy a F = &b ]n[2S(t — 1) + 1] and an AF segment 
a A = {£ odd) and fc B ln(2S + 1) (£ even). The total 
entropy contribution per site is 



Act = 



N 



(3.8) 



where the sum is over all the segments and Oi is the en- 
tropy of segment i as discussed above. We note that 



N = J2i n i- Splitting the sums into sums over AF seg- 
ments and FM segments and using the fact that the num- 
ber of AF segments equals the number of FM segments, 
we obtain 



Act 



(a a) + (gg) 
(n A ) + (n F ) ' 



(3.9) 



where (<j a ,f) is the average entropy of an AF or FM seg- 
ment and (tia,f) is the average number of sites in a seg- 
ment. For reasons discussed above the appropriate way 
to count boundary sites between FM and AF segments is 
to attribute them to the AF side. Therefore n F =1 — 1 
and ha = & + 1- The probability to find a FM or AF 
segment with £ bonds is 



P F (£) = (1 - v ) v l - x (FM segments), 
P A {£) = p(l - p)^ 1 (AF segments), 



(3.10) 



(3.11) 



respectively, and we can calculate the desired averages as 

oo 

{o f ) = Y j Pf{() ln[25(£-l) + l], 

i=i 

oo 

(CT A )=ln(25+l)^P A (2m), 

m— 1 

oo 

(n F ) = J2MW-l) = T J ^, 

OO 

(n A ) = ^PAme +!) = —■ 

In Fig. ^ we plot Act as a function of p, which agrees 
quite well with the estimates of the missing entropy from 
the HTE data. 

Let us now turn to the Curie behavior of the suscepti- 
bility in the intermediate-temperature regime. The idea 
of nearly independent effective spin degrees of freedom 
implies, of course, a 1/T-dependence of the susceptibil- 
ity as long as they are uncorrelated. Analogous to the 
entropy we can calculate the effective Curie constant c 
by averaging the effective spin sizes (x = p . c/fcsT): 



1 (S 2 F ) + (S 2 A ) 
3 (n F ) + (n A ) 



(3.12) 



with 



(S F ) = J2P F (£)S(£-1)[S(£-1) + 1], 

1=1 

OO 

m— 1 

which leads to 

_ P S[S(2p 3 - Ap 2 + p - 1) +p - 1] 
C ~ 3(l-p)(p-2) • 



(3.13) 



(3.14) 



G 



At p = 0.5 we obtain c = | for S — ^ and c = ^ 
(= 0.3888) for S = 1, both of which agree well with the 
values deduced from the HTE and the TMM, c = 0.13 
(S = |) and c= 0.4 (5 = 1). 

This discussion of our data is clearly consistent with 
the picture of weakly coupled FM and AF segments. The 
intermediate temperature regime is aparently well de- 
scribed by an ensemble of independent AF and FM finite 
size spin chains with a certain probability distribution. 



IV. LOW-ENERGY EXCITATIONS 

In the previous section we based the interpretation of 
our data on the argument that effective spin degrees of 
freedom are formed on FM and AF segments which cou- 
ple only weakly among each other. In this section we 
would like to substantiate this picture by discussing some 
typical cases and using data obtain from exact diagonal- 
ization. For simplicity we will concentrate on the spin-^ 
system. 

Let us first consider the case of small p where almost 
all bonds are antiferromagnetic. This system consists of 
very long AF segments separated by mostly single FM 
bonds. The extreme case of one FM bond between two 
semi-infinite AF segments corresponds to a FM impurity 
bond in an AF spin-1 chain. By using the bosoniza- 
tion and the RG techniques, a weak FM coupling can 
be shown to.-he an irrelevant perturbation with scaling 
dimension 2 £3 Although the FM coupling of strength J 
in our model is not weak, it is very likely that, it will 
be renormalized to zero in the low-energy limit EB Hence 
the FM bond surrounded by infinitely long AF segments 
will completely decouple the AF segments. For finite but 
small p the finite lengths of the AF segments stops the 
renormalization of the FM bonds at some energy scale 
determined by I ~ 1/p. For p <C 1 the renormalized 
FM couplings are so small that we can neglect them in 
the first approximation. A decoupled AF segment of even 
length form a singlet ground state while a segment of odd 
length form a doublet with spin i. Thus, in this approxi- 
mation the degenerate ground state consists of decoupled 
AF segments of even and odd length carrying spin and 
1, respectively. The residual weak FM couplings will 
introduce small nearest-neighbor couplings between the 
doublets and lift the degeneracy. Since the even-length 
segments forms singlets, they do not contribute to the 
low-energy degrees of freedom, but they are important to 
determine the effective couplings between the spin-^ seg- 
ments; two spin-i segments separated by an even num- 
ber (including 0) of singlets will couple ferromagnetically 
while an odd number of separating singlets makes the ef- 
fective coupling antiferromagnetic. The actual value of 
the effective coupling depends on the length of the two 
spin-i segments as well as the number of and lengths of 
the separating singlets. Therefore the resulting effective 
Hamiltonian is again that of a random spin-i Heisenberg 



chain but with random bond strength of either sign. 

We have confirmed these properties by exact diago- 
nalization of finite spin chains for some typical segment 
configurations. Figure || shows the segmentation of a 
chain with two AF segments of odd lengths 5 and 7 sep- 
arated by a single FM bond. In the spectrum the two 
lowest energy states, a singlet and a triplet, are very 
close in energy (AE$ = 0.13J), while the gap between 
the ground state and the next excited state which in- 
volves intra-segment excitations is considerably higher, 
AEm = 0.56J. This demonstrates the separation of the 
low-energy scale from the intermediate one. The size of 
the splitting between the lowest singlet and triplet and 
the relative location of these two levels define the magni- 
tude and sign of the effective coupling: J cff = —0.13 J. In 
Fig. ^| we show a chain with two segments of odd num- 
ber of bonds separated by one singlet and the data of 
their energy levels. In this case the corresponding en- 
ergy gaps are AE S — 0.047J and AE M = 0.65J, and 
the separation of energy scales is even more pronounced. 
The relative position of the low-energy singlet and triplet 
confirms that the effective coupling in this case is anti- 
ferromagnetic. In Fig. [l^ we link the two configurations 
in Figs. |^ and || together. The spectrum of the effective 
Hamiltonian, 

H cS = E a + jfSi • S 2 + J 2 cff S 2 • S 3 (4.1) 

with J? s = 0.047 J and Jf = -0.125 J estimated from 
Figs, g and ^| (dashed lines), agrees very well with the 
lowest part of the exact spectrum (solid lines), show- 
ing that the low-energy degrees of freedom are well de- 
scribed by a nearest-neighbor Heisenberg Hamiltonian. 
The energy shift Eq = — 7.65J is adjusted so that the 
ground-state energies agree. The fact that the true spec- 
trum is reproduced so well with effective coupling cal- 
culated from two-segment chains confirms that possible 
non-nearest-neighbor interactions in the effective Hamil- 
tonian are small. 

Let us now turn to the case of p close to 1 , the nearly 
FM chain. Although different in structure, we can argue 
quite similarly to p C 1 on the separation of the low- 
energy scale. Consider two FM segments of finite length 
separated by a single AF bond. This AF bond has a ten- 
dency to lock its two adjacent spins into a singlet, which 
effectively decouples the two FM segments. In a mod- 
ified chain where all AF bonds are much stronger than 
the FM ones, this tendency is even more pronounced, 
and in the limit of infinite AF exchange the FM seg- 
ments are completely decoupled. In this limit each FM 
segment forms a local ground state of maximum spin, 
and the overall ground state is highly degenerate since 
the large segment-spins are non-interacting. A finite Jaf 
lifts this degeneracy, and the effective Hamiltonian for 
the segment-spins is, to first order in Jfm/Jaf, an AF 
nearest-neighbor Heisenberg Hamiltonian with coupling 
strengths depending on the lengths of the FM segments. 
Our exact diagonalization results suggest that this pic- 
ture remains true even when Jfm/Jaf = 1 and that 
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this results in a separation of energy scales. In Fig. [If] 
is a FM chain with a single AF bond. From the seg- 
mentation we expect five almost degenerate multiplets 
of spin S — i , I , . . . , I corresponding to different rela- 
tive orientations between the effective spins Si = 2 and 
§2 = I- These multiplets indeed appear in the spectrum 
[Fig. [ll](b)] and are reproduced by the effective Hamilto- 
nian 



H = E a 



(4.2) 



with E = -2.72J, J x = 8.5 x 1CT 3 J, J 2 = 3.2 x 1(T 4 J, 
J 3 = 2.0 x 10~ 5 J, and J 4 = -7.8 x 1(T 7 J. In the spec- 
trum in Fig. [ll](b) we can also clearly identify the four 
multiplets S = \ , . . . , | corresponding to the right seg- 
ment forming a spin-| state (magnon excitation) and the 



three multiplets S 



3 5 



and | corresponding to the 



left segment forming a spin-1 state. The excitation ener- 
gies to these one-magnon states are AEj^ = 0.22J and 
AEm 2 = 0.31J, respectively. Finally we form a three- 
segment chain in Fig. The lower part of the true 
spectrum is once again in excellent agreement with the 
spectrum of the effective Hamiltonian which is taken to 
be 



T~Lcfi — Eq + Ji ff Si • S2 ■ 



J 2 off S 2 



•So 



(4.3) 



In Eq. Q Jf = 0.0085J is estimated from Fig. [n] 
and Jf ff = 0.011J is obtained from a similar calculation. 
The energy shift is adjusted to Eq = — 4.28J for a best 
fit. These results confirm that the low-energy physics 
for spin chains with p close to one is also well described 
by a random nearest-neighbor Heisenberg Hamiltonian. 
However, in contrast to small p the effective couplings 
for p close to one arc all antifcrromagnctic, though the 
magnitude of the spins are random. 

For intermediate values of p the situation is more com- 
plex because the typical configuration of FM and AF 
bonds leads to a sequence of rather short segments. The 
extreme case is a sequence of alternating FM and AF 
bonds which is known to favor a singlet (dimer) configu- 
ration with an excitation gap. Such a sequence, located 
between two segments with finite effective spin degrees of 
freedom, yields extremely small couplings. On the other 
hand, we can also find the situation that a FM segment 
borders on an AF segment with an even number of bonds, 
which forms an effective spin doublet. In this case the 
effective coupling between the two segments is not neces- 
sarily small. Nevertheless, the coupling is small enough 
to form the separate segment spins which in turn gov- 
ern the low-energy physics. This type of configurations 
would, of course, be the first to undergo intersegment 
correlation as the temperature is lowered. Hence, we ex- 
pect that in the effective low-energy Heisenberg model 
the distribution of the coupling strengths as well as spin 
sizes is rather broad for these systems. 



V. CONCLUSIONS 



We have studied the ID quantum Heisenberg model 
with random FM and AF bond disorder by means of 
several numerical methods. A comparison of the quan- 
tum spin chains with the corresponding classical system 
reveals interesting similarities. At first sight the differ- 
ences seem to be only of quantitative nature. However, 
through a careful analysis of our data, we found that 
the low-temperature (low-energy) physics of the quantum 
and classical system has profound differences. While the 
quantum spin system exhibits a separation of two energy 
scales between the low and intermediate energies, only 
one scale is present in the classical case. 

The idea of the creation of weakly-interacting effective 
spin degrees of freedom in FM and AF segments along 
the chain plays a key role in understanding the physics 
of the quantum system. While the high-energy (high- 
temperature) behavior is determined by the individual 
spins of the chain, they are replaced by new effective 
spin degrees of freedom towards the low-energy limit. As 
we have demonstrated in the last section, the low-energy 
physics is well captured by an effective model which has 
the structure of a nearest-neighbor Heisenberg chain with 
spins of variable size and FM and AF interactions of 
random strength. As a consequence three temperature 
regimes emerge in our quantum spin system (in contrast 
to two in the classical case): the high-, intermediate- and 
low-temperature regime. The latter two are clearly de- 
termined by the effective low-energy model. In this sense 
the intermediate temperature regime may be considered 
as the high-temperature regime of the effective model. 
Therefore the effective spins yield a Curie-like suscepti- 
bility in this regime as confirmed in Section IIIC. In a 
similar way the specific heat may be separated into two 
contributions. One is due to the individual spins, which 
generate a specific heat peak when they start to correlate 
at a temperature fc^T ~ J. The other originates from 
the effective spins, which give rise to another presumably 
much weaker peak-like structure containing, however, a 
considerable fraction of the entropy at lower temperature. 

The methods we used in this study did not allow us to 
investigate the low-energy physics extensively. Alterna- 
tive methods have to be applied here. The investigation 
of the thermodynamics and magnetic properties of this 
regime is in progress and will be presented elsewhere. 
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FIG. 1. Specific heat of the spin-i chain for p = (solid), 
0.25 (dotted), 0.5 (long dashed), 0.75 (short dashed), and 1 
(dashed dotted) obtained from two-point Pade approxima- 
tion. The filled circles denote the results by Blote for the 
uniform chains and the empty circles the data of our transfer 
matrix calculation. 



FIG. 2. Specific heat of the spin-1 chain for p — (solid), 
0.25 (dotted), 0.5 (long dashed), 0.75 (short dashed), and 
1 (dashed dotted) obtained from Dlog Pade approximation. 
The filled circles denote the results by Blote for the uniform 
chains and the empty circles the data of our transfer matrix 
calculation. 
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FIG. 3. Specific heat divided by temperature for the spin-i 
chain. The notations are the same as in Fig. [j]. 



FIG. 4. Specific heat divided by temperature for the spin-1 
chain. The notations are the same as in Fig. ^. 



FIG. 5. Inverse susceptibility of the spin-^ chain for p — 0, 
0.25, 0.5, 0.75 and 1. 



FIG. 6. Inverse susceptibility of the spin-i and spin-1 
chains at p = 0.5. 



FIG. 7. Missing entropy. The solid (S = \) and d ashed 



(S — 1) curves are calculated from Eqs. (3.9) and (3.11). The 
circles (diamonds) are the missing entropies estimated from 
the HTE for S= \ (S = 1). 



FIG. 8. (a) The segmentation of a short chain containing 
one FM bond. AF couplings are represented with + and FM 
couplings with — . (b) The energy spectrum of the chain in 
(a). 



FIG. 9. (a) The segmentation of a chain with two FM 
bonds resulting in two odd-length segments separated by one 
singlet. Bonds are represented as in Fig. ^. (b) The energy 
spectrum of the chain in (a). 



FIG. 10. (a) The segmentation of a chain of length L = 19. 
The segments appearing in this particular chain are the ones 
in Figs. ^ and || Bonds are represented as in Fig. ^. (b) 
The energy spectrum of the chain in (a) (soild lines) and the 
spectrum of the correspinding effective Hamiltonian for the 
segment-spins (dashed lines). 



FIG. 11. (a) An AF bond in a FM chain locks its two ad- 
jecent spins into a singlet. This effectively decouples the FM 
segments which form local ground states of largest possible 
spins. Bonds are represented as in Fig. ||. (b) The spectrum 
of the chain in (a). 



FIG. 12. (a) The segmentation of a nearly FM chain result- 
ing in three weakly interacting large segment-spins. Bonds are 
represented as in Fig. ^. (b) The energy spectrum of the chain 
in (a) (solid lines) and the spectrum of the corresponding ef- 
fective Hamiltonian (dashed lines). 
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TABLE I. 


Coefficients u(l,m 


) for the spin-i chain; u(K;p) 


= ( J /4) Ei E 


u(l,m){pJ/4) l p m . 


(1, m) 


u(l, m) 


(I, m) 


u(l, m) 


(/, m) 


u(l, m) 


(1,0) 


-3.0 


(12,3) 


7713.500413 


(18,3) 


-1134121.9573 


(2,0) 


-3.0 


(13,0) 


8196.990854 


(18,4) 


-288649.3274 


(2,1) 


6.0 


(13,1) 


-25934.104588 


(18,5) 


115459.7310 


(3,0) 


5.0 


(13,2) 


23117.786328 


(19,0) 


-2079.9826 


(4,0) 


15.0 


(13,3) 


5632.636519 


(19,1) 


-646413.7595 


(4,1) 


-30.0 


(13,4) 


-2816.318258 


(19,2) 


2645440.1027 


(5,0) 


-4.2 


(14,0) 


6660.909634 


(19,3) 


-3976953.5968 


(5,1) 


-20.8 


(14,1) 


-39516.948771 


(19,4) 


1935729.0738 


(5,2) 


20.8 


(14,2) 


78429.276975 


(19,5) 


63297.2706 


(6,0) 


-61.13333333 


(14,3) 


-51765.812859 


(19,6) 


-21099.0874 


(6,1) 


122.26666667 


(14,4) 


-780.557689 


(20,0) 


-1251808.4999 


(7,0) 


-40.48571429 


(14,5) 


312.223089 


(20,1) 


5117166.1965 


(7,1) 


201.75238095 


(15,0) 


-28013.413242 


(20,2) 


-7236338.9668 


(7,2) 


-201.75238095 


(15,1) 


80562.547491 


(20,3) 


2810586.2334 


(8,0) 


204.90476190 


(15,2) 


-47628.744538 


(20,4) 


3017931.4549 


(8,1) 


-381.46666667 


(15,3) 


-65867.605907 


(20,5) 


-1202116.2591 


(8,2) 


-85.02857143 


(15,4) 


32933.802930 


(20,6) 


-5056.3158 


(8,3) 


56.68571428 


(16,0) 


-54542.64412 


(20,7) 


1444.6371 


(9,0) 


353.8973545 


(16,1) 


255075.95809 


(21,0) 


-1184631.917 


(9,1) 


-1276.5291005 


(16,2) 


-433951.49101 


(21,1) 


7782735.876 


(9,2) 


1276.5291005 


(16,3) 


275899.26561 


(21,2) 


-19645416.929 


(10,0) 


-500.5733333 


(16,4) 


20102.59259 


(21,3) 


23366818.910 


(10,1) 


571.1106878 


(16,5) 


-8041.03709 


(21,4) 


-10787051.462 


(10,2) 


1290.1079365 


(17,0) 


66386.55357 


(21,5) 


-1075629.595 


(10,3) 


-860.0719576 


(17,1) 


-123625.66748 


(21,6) 


358543.182 


(11,0) 


-1918.0408850 


(17,2) 


-161124.99951 


(22,0) 


4217217.366 


(H.1) 


6375.6162771 


(17,3) 


568820.19930 


(22,1) 


-15639609.441 


(H,2) 


-6238.2222222 


(17,4) 


-282707.26277 


(22,2) 


16514166.308 


(11,3) 


-274.7881097 


(17,5) 


-2043.40448 


(22,3) 


5968948.429 


(H,4) 


137.3940551 


(17,6) 


681.13439 


(22,4) 


-25376121.949 


(12,0) 


259.304916 


(18,0) 


293845.3269 


(22,5) 


9967514.773 


(12,1) 


3338.140375 


(18,1) 


-1270211.3634 


(22,6) 


182933.980 


(12,2) 


-11570.250620 


(18,2) 


1989832.2634 


(22 7) 


-52266.733 
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TABLE II. Coefficients x(/,m) for the spin-± chain; x{K;p) = /3fi 2 J^i E m X(l,m)(/3J / '4)9 



(Z, m) 




(l,m) 


X{l,m) 


(l,m) 


X{l,m) 


(0,0) 


0.25 


(6,6) 


32.0 


(9,9) 


256.0 


(1,0) 


-0.5 


(7,0) 


1.015873016 


(10,0) 


-59.9102293 


(1,1) 


1.0 


(7,1) 


8.679365080 


(10,1) 


24.7207055 


(2,1) 


-2.0 


(7,2) 


-0.133333333 


(10,2) 


-134.8949559 


(2,2) 


2.0 


(7,3) 


-74.577777778 


(10,3) 


-241.7721340 


(3,0) 


0.666666667 


(7,5) 


224.0 


(10,4) 


508.1876543 


(3,1) 


0.666666667 


(7,6) 


-224.0 


(10,5) 


1173.6380953 


(3,2) 


-6.0 


(7,7) 


64.0 


(10,6) 


-1585.8793651 


(3,3) 


4.0 


(8,0) 


18.12857143 


(10,7) 


-2048.0 


(4,0) 


0.833333333 


(8,1) 


-7.32063492 


(10,8) 


4352.0 


(4,1) 


2.0 


(8,2) 


57.50476191 


(10,9) 


-2560.0 


(4,2) 


6.0 


(8,3) 


38.29841270 


(10,10) 


512.0 


(4,3) 


-16.0 


(8,4) 


-237.81587302 


(n,o) 


-109.9150746 


(4,4) 


8.0 


(8,5) 


-96.0 


(11,1) 


300.2769921 


(5,0) 


-1.4 


(8,6) 


629.33333333 


(H,2) 


70.4828219 


(5,1) 


-1.866666667 


(8,7) 


-512.0 


(11,3) 


-338.9203527 


(5,2) 


6.0 


(8,8) 


128.0 


(H,4) 


-1239.6134039 


(5,3) 


22.666666667 


(9,0) 


13.1816578 


(11,5) 


983.6675837 


(5,4) 


-40.0 


(9,1) 


-52.9643739 


(11,6) 


4324.9777778 


(5,5) 


16.0 


(9,2) 


-23.9428571 


(11,7) 


-3283.7079367 


(6,0) 


-4.433333333 


(9,3) 


180.4486772 


(11,8) 


-6912.0 


(6,1) 


-0.888888889 


(9,4) 


259.9365079 


(11,9) 


10922.6666667 


(6,2) 


-20.444444444 


(9,5) 


-656.5079365 


(11,10) 


-5632.0 


(6,3) 


10.666666667 


(9,6) 


-522.6666667 


(11,11) 


1024.0 


(6,4) 


74.666666667 


(9,7) 


1685.3333333 






(6,5) 


—96.0 


(9,8) 


-1152.0 






^ V^ RT TTT 
±l\DljFj 111. 


Coefficients u{1,tti 


) for the spin-1 chain; u(K;p) — 






(l,m) 


u(l, m) 


(I, m) 


u(l, m) 


(I, m) 


u(l, in) 


(1,0) 


-2.6666666667 


(10,2) 


28.8792214 


(15,1) 


45953.430459 


(2,0) 


-1.3333333333 


(10,3) 


-19.2528142 


(15,2) 


-45757.817154 


(2,1) 


2.6666666667 


(11,0) 


481.3638086 


(15,3) 


-391.226630 


(3,0) 


5.9259259259 


(11,1) 


1301.8758468 


(15,4) 


195.613379 


(4,0) 


6.6666666667 


(11,2) 


-1301.3134946 


(16,0) 


54519.198734 


(4,1) 


-13.333333333 


(11,3) 


-1.1247047 


(16,1) 


-99299.036124 


(5,0) 


-18.251851852 


(H,4) 


0.5623535 


(16,2) 


-29214.214665 


(5,1) 


-1.896296296 


(12,0) 


2977.0623249 


(16,3) 


19463.245209 


(5,2) 


1.896296296 


(12,1) 


-5835.6140989 


(16,4) 


19.346878 


(6,0) 


-31.664197531 


(12,2) 


-355.5316528 


(16,5) 


-7.738804 


(6,1) 


63.328395062 


(12,3) 


237.0211016 


(17,0) 


28441.68270 


(7,0) 


58.441544189 


(13,0) 


-849.9549193 


(17,1) 


-249359.16806 


(7,1) 


23.237154614 


(13,1) 


-8019.3116116 


(17,2) 


246954.37836 


(7,2) 


-23.237154614 


(13,2) 


8006.2448768 


(17,3) 


4809.45489 


(8,0) 


147.08759553 


(13,3) 


26.1334724 


(17,4) 


-2404.41627 


(8,1) 


-293.62210464 


(13,4) 


-13.0667454 


(17,5) 


-0.37317 


(8,2) 


-1.65925926 


(14,0) 


-12916.439014 


(17,6) 


0.12400 


(8,3) 


1.10617283 


(14,1) 


24667.952605 


(18,0) 


1688755.7953 


(9,0) 


-179.33842257 


(14,2) 


3494.640784 


(18,1) 


-5177380.6581 


(9,1) 


-190.19539704 


(14,3) 


-2329.308912 


(18,2) 


5503050.1273 


(9,2) 


190.19539704 


(14,4) 


-0.677419 


(18,3) 


-4013509.8333 


(10,0) 


-669.2886525 


(14,5) 


0.270973 


(18,4) 


517214.6225 


(10,1) 


1328.9508980 


(15,0) 


-1530.190525 


(18,5) 


-206885.8485 
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TABLE IV. Coefficients x(J,m) for the spin-1 chain; x(K;p) = /V J^i J2 m m )(/ 3J / 2 )V' 



(I, m) 


X(l,m) 


(l,m) 


X(l,m) 


(l,m) 




(0,0) 


0.6666666667 


(5,3) 


299.19341564 


(8,0) 


-5.234080 


(1,0) 


-1.7777777778 


(5,4) 


-449.49245542 


(8,1) 


6.320244 


(1,1) 


3.5555555556 


(5,5) 


179.79698217 


(8,2) 


425.727623 


(2,0) 


1.4814814815 


(6,0) 


4.14229538 


(8,3) 


1053.738742 


(2,1) 


-9.4814814815 


(6,1) 


-8.95473251 


(8,4) 


-1911.972048 


(2,2) 


9.4814814815 


(6,2) 


-148.3676269 


(8,5) 


-7884.430625 


(3,0) 


1.580246914 


(6,3) 


-164.8139003 


(8,6) 


18539.066606 


(3,1) 


9.481481481 


(6,4) 


1281.0534979 


(8,7) 


-13637.934055 


(3,2) 


-37.925925926 


(6,5) 


-1438.3758573 


(8,8) 


3409.483514 


(3,3) 


-25.283950617 


(6,6) 


479.4586191 


(9,0) 


-67.132623 


(4,0) 


-2.13991770 


(7,0) 


21.0825005 


(9,1) 


-461.172347 


(4,1) 


6.32098765 


(7,1) 


50.8237507 


(9,2) 


2085.379854 


(4,2) 


61.10288066 


(7,2) 


85.6220698 


(9,3) 


— 10591.213175 


(4,3) 


-134.84773663 


(7,3) 


-633.0976376 


(9,4) 


42515.683875 


(4,4) 


67.42386831 


(7,4) 


-1373.4491693 


(9,5) 


-74434.761484 


(5,0) 


-5.37283951 


(7,5) 


5024.3267795 


(9,6) 


57428.487934 


(5,1) 


-19.45459534 


(7,6) 


-4474.9471117 


(9,7) 


-16408.139410 


(5,2) 


0.70233196 


(7,7) 


1278.5563176 







TABLE V. The ground state energy estimated from Dlog Pade approximants. 



p 





0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 


0.9 


1 


5=1/2 
S = 1 


-0.439 
-1.395 a 


-0.428 
-1.35 


-0.414 
-1.308 


-0.398 
-1.267 


-0.380 
-1.222 


-0.360 
-1.183 


-0.337 
-1.143 


-0.316 
-1.104 


-0.293 
-1.06 


-0.273 
-1.04 


-0.240 
-1.01 



a This value is in reasonable agreement with the result obtained by Nightingale and Blote: E3 E a = -1.4015 J. 
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